Random walk and Pair- Annihilation Processes on Scale-Free Networks 
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We investigate the dynamic scaling properties of stochastic particle systems on a non-deterministic 
scale-free network. It has been known that the dynamic scaling behavior depends on the degree 
distribution exponent of the underlying scale-free network. Our study shows that it also depends 
on the global structure of the underlying network. In random walks on the tree structure scale-free 
network, we find that the relaxation time follows a power-law scaling r ~ N with the network size 
N. And the random walker return probability decays algebraically with the decay exponent which 
varies from node to node. On the other hand, in random walks on the looped scale-free network, 
they do not show the power-law scaling. We also study a pair-annihilation process on the scale-free 
network with the tree and the looped structure, respectively. We find that the particle density 
decays algebraically in time both cases, but with the different exponent. 



I. INTRODUCTION 

Networks composed of nodes and edges have been at- 
tracting a lot of research interest recently Q. Unlike 
conventional networks such as a periodic regular network 
and a random network, many real-world networks have 
complex structure. Many of them belong to the class of 
so-called scale-free (SF) networks. Denoting the degree 
k of a node as the number of edges attached to it, a SF 
is characterized with the power-law degree distribution 



(1) 



Here 7 is called the degree distribution exponent. 

The power-law degree distribution implies that the SF 
network has an inhomogeneous structure. On the one 
hand, it is a challenging problem to characterize and un- 
derstand the organization principle of the SF network. 
On the other hand, it is also interesting to study thermo- 
dynamic or dynamic systems on such an inhomogeneous 
structure. Various physical problems have been studied. 
Examples include the ferromagnetic phase transitions in 
the Ising model [2( , the non-equilibrium phase transition 
in the epidemic spreading model the random walk 
process H H H HJi, the dynamic scaling m pair- 
annihilation process |10t fllj . and so on. In those studies, 
the 7-dependent scaling properties have been studied. 

In the present work, we address the question how the 
global structure of the underlying SF network influences 
the dynamic scaling behavior of stochastic particle sys- 
tems. For that purpose we study the random walk pro- 
cess and the pair-annihilation process, and investigate 
their dynamic scaling property. Both systems are studied 
on the tree structure SF (TSF) networks and the looped 
structure SF (LSF) networks, respectively. 

To be specific, we study the stochastic systems on 
the Dorogovtsev-Mendes-Samukhin (DMS) network [l2| . 
which is a generalization of the Barabasi-Albert (BA) 
network It is a non- deterministic model for a grow- 
ing network: Each time step, a new node is added, and 
linked with m nodes which are selected among existing 
nodes with the probability given by ITj oc (fcj + a). Here 



the parameter a is called an initial attractiveness. The 
BA model corresponds to the a = case of the DMS net- 
work. The resulting network is scale-free and the degree 
distribution exponent is given by 
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(2) 



With the parameters m and a, one can vary the value 
of the degree exponent. At the same time, one can also 
generate a TSF network with m = 1 or a LSF network 
with ro^l. Hence, we can study the effect of the degree 
distribution exponent and the global network structure 
systematically. 

This paper is organized as follows: In Sec. [HJ we 
present the results on the random walk process. Sec- 
tion IIIII is devoted to the scaling property of the pair- 
annihilation process. Summary will be given in Sec. IIVI 



II. RANDOM WALK PROCESS 

As a basic and fundamental stochastic process, the ran- 
dom walk process on complex networks has been at- 
tracting a lot of research interest 

11 nasi in 

this work we concentrate on relaxation dynamics on SF 
networks. 

On a network with N nodes, a random walker is as- 
signed to a starting node denoted by s at time t = 0. 
Then, at each unit time step At = 1, it hops to one of 
the neighboring nodes selected randomly with the equal 
probability. Defining P(i, t; s, 0) as the probability to find 
the walker at node i at time t, one finds that it evolves 
in time as 



P(i,t+l; S ,0) = ^4^ P (^ s <°) 



3=1 



(3) 



with the initial condition P(z,0;s,0) = <5j iS . Here A = 
{Ay} is the adjacency matrix whose elements are Aij = 
1 (0) if two nodes i and j are connected (disconnected) . 

The relaxation dynamics is studied with the so-called 
return probability R s (t) = P(s,i;s,0). For a given net- 
work, one can solve numerically the master equation 




10 10 10 10 10 10 10 




OON = 


1000 


cxiN = 


4000 




16000 


e^N = 


64000 



FIG. 1: The return probability Rh (open symbols) and 
Rp (rilled symbols) in the DMS networks with 7 = 3 and 
with m = 1 in (a) and m — 2 in (b). 




FIG. 2: Relaxation time r for Rij(t) (open symbols) and 
Rp(t) (filled symbols) on the tree structure SF networks. The 
solid line has the slope 1. 



Eq. J2J iteratively to obtain R s (t). We then average it 
over different realizations of the networks. 

On a network with loops, the probability distribution 
converges to the stationary one P s tat.{i) — fei/QDj ^i) m 
the t -> oo limit [(|. Hence, the return probability con- 
verges to R s (t = oo) = k s /(J2j kj). On the other hand, 
on a tree network, the probability distribution is oscil- 
lating in time, so is the return probability. The random 
walker cannot return to a starting node in odd time steps, 
which means that R s (t) = at odd t. In this case, we 
only measure the return probability at even time steps, 
which converges to stationary value 2k s /(J2j kj). 

In a SF network, the node degree is distributed so 
broadly that one can define an exponent q s for each node 
s describing the degree scaling with the network size: 



(4) 



For instance, in the DMS network with the degree distri- 
bution exponent 7, a peripheral node with the minimum 
degree has q = 0, while the hub with the maximum de- 
gree has q = 1/(7— 1). Then, the return probability in 
the stationary state scales with the network size as 



R s {t^ 00) ~ aMi-9=) . 



(5) 



Due to the broad degree distribution, the return prob- 
ability R s (t) may have a different scaling property at 
different starting node s with different values of q s . So, 
we measured the return probability Rn(t) for the hub 
and Rp(t) for a peripheral node on the DMS networks. 
We present the numerical data for the return probability 
Rh and Rp in Fig. ^ We compare the data obtained 
from the TSF networks and the LSF networks. The data 
show that the return probability behaves distinctly de- 
pending on the global structure of the network, which 
will be discussed in detail in the following. 



A. Random walks on tree networks 

We perform the quantitative analysis of the scaling 
property in the TSF. From Fig. ^ one finds that the 



return probability relaxes much slower in the TSF net- 
work with m — 1. We estimated the relaxation time 
r using the condition R s (t = r) = cR s (t = 00) with 
a constant c — 2. The relaxation time is plotted as a 
function of the network size N in Fig. El It shows that 
the relaxation time estimated from Rh and Rp follows 
a power-law scaling 



r ~ N z 

with the same dynamic exponent 
z = 1.0 . 



(6) 



(7) 



We provide a theoretical argument for the numeri- 
cal result from the analysis of the mean first passage 
time (MFPT). The MFPT problem has been studied on 
complex networks recently 0, 0, IUj an( i some rigorous 
results are known [HQ- Consider the MFPT, denoted 
by Tj t i, from an arbitrary node i (degree fcj) to one of 
its neighboring node j. Let {ni, n%, ■ • ■ , = j} denote 
the neighbors of the node i. Then, following Ref. 0, the 
MFPT on tree networks satisfies the recursion relation 



T, 



5> 



(8) 



Applying the recursion relation repeatedly until one ar- 
rives at dangling ends, one can find the explicit solution 
for Tj^. Without the link between i and j, the tree net- 
work would be decomposed into two parts. Denoting the 
number of nodes in the i side by Ni, the MFPT is simply 
given by 



(9) 



Hence, for a typical adjacent nodes i and j, one has that 
Tj t i ~ iV 1 . For non-adjacent nodes i and j, the MFPT is 
given by the sum of the MFPT's given by Eq. 10 along 
the path between them. Hence, the MFPT between a 
typical node pair is given by T ~ D^N, where Z?at is 
the mean diameter of the networks. For SF networks, 
the mean diameter scales at most logarithmically with 



3 




tIN r/N 



FIG. 3: Scaling analysis of Rh (open symbols) and Rp (filled 
symbols) in the TSF networks with 7 = 3 (a) and 7 = 4 (b). 



the network size |l5j. Therefore, we conclude that the 
relaxation time follows the power-law scaling r ~ N z 
with the dynamic exponent z = 1. 

The power-law scaling of the relaxation time suggests 
that the return probability decays algebraically as 

R s (t)~t- S ° (10) 

with the decay exponent S s . For a node s with the degree 
scaling k s ~ N qs , the decay exponent can be deduced 
from the finite-size-scaling ansatz 

R s (t) = N-^-^f(t/N z ) (11) 

with the dynamic exponent 2 = 1. For large t N z , 
the return probability should converge to the stationary 
value R s (t = 00) ~ N-^-i^) ( S ee Eq. (J5J). It requires 
that the scaling function should behave as f(x ^> 1) = 
constant. The power-law scaling for t <C N z requires that 
the scaling function should behave as f(x < 1) ~ x~ ds , 
which yields that R s ~ t -S, N s e z-(i-q s ) _ Therefore, the 
finite-size-scaling ansatz predicts that the decay expo- 
nent is given by 

S s = l-q s . (12) 

It is interesting to note that the decay exponent varies 
with the degree scaling exponent q a of the starting node 
s. We confirm the scaling behavior with the scaling plot 
of i? s iV 1_9s versus t/N for the hub and the peripheral 
node in Fig. E| For 7 = 3, q H = 1/(7-1) = 1/2 for 
the hub and qp = for the peripheral node, which yields 
that 8h = 1/2 and Sp = 1. Similarly, for 7 = 4, one 
expects that 5h = 2/3 and 5p = 1. One finds that all 
data from different network sizes collapse very well in the 
scaling plot, which supports the result in Eq. (|T2*)l . 

It is also interesting that the return probability decays 
as R s (t) ~ t~ Sa with the exponent S s < 1 at all nodes. 
It indicates that the random walks are recurrent in the 
N -> 00 limit [l4| . 




l IN 1 / In N 

FIG. 4: Plots of 6 obtained from R P (t) on the DMS networks 
with 7 = 2.5, 3, 4 and with m — 2 (open symbols) and m = 
4 (filled symbols). # is plotted against 1/N in (a) and 1/lniV 
in(b). 

B. Random walks on looped networks 

From Fig. ^ one can see that the return probability 
decays much faster in the LSF networks. The downward 
curvature in the log-log plot implies that the decay is 
faster than a power-law decay. The relaxation time mea- 
surement also indicates a faster decay. We found that 
the relaxation time, measured using the condition that 
R s (t) = 2i? s (oo), grows at most logarithmically with 
the network size N. This is in contrast to the power- 
law growth in the TSF networks. In this subsection, we 
address a question how the return probability decays in 
time in looped SF networks. 

In the random network, it is known that the return 
probability follows a stretched exponential decay as 

R s {t) - R s (oo) ~ e~ ate (13) 

with a constant a and the exponent = 1/3 |l6j . It was 
reported that the return probability in the small-world 
network follows the stretched exponential decay, too ^3 • 
It might suggest that the looped SF network follow the 
stretched exponential decay, too. 

On the DMS networks with m > 2, we measured R s (t) 
numerically and fitted the data to the form in Eq. (|13|) 
to estimate 8. If the return probability follows an expo- 
nential decay, one would have 6 = 1. On the other hand, 
one would have 9 < 1, if it follows the stretched expo- 
nential decay. In Fig.QJ we present the data for 0(N) for 
Rp(t) obtained on the DMS networks with various values 
of 7 = 2.5, 3, 4 and m = 2, 4 up to sizes N ~ 10 6 . Similar 
behaviors are observed in 9(N) for Rn(t). At small val- 
ues of N, the exponent seems to depend on m and to be 
smaller than 1. One may interpret that as an evidence 
of the stretched exponential decay with a non-universal 
exponent 9. 

However, we observe that there is a very strong finite 
size effect. The plot of 9 versus 1/ In in Fig.0|(b) shows 
that the finite size effect is significant even for N ~ 10 6 . 
With this strong finite size effect, one can not exclude the 
possibility of the exponential decay with 9 = 1. We sug- 
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gest that an analytic approach be necessary to conclude 
whether the return probability follows the exponential or 
the stretched exponential decay. 

Before closing this section, we remark on the previ- 
ous studies on the random walks on deterministic SF 
networks 0, ■ Unlike the non-deterministic LSF net- 
works studied in this work, the deterministic LSF net- 
works [la, [l9|, |20( behave similarly as the TSF networks. 
For example, in the hierarchical network 20] which has a 
looped structure, the relaxation time scales algebraically 
as r ~ N 0. The reason why the hierarchical network 
behaves as TSF networks is clear. The network has a 
symmetry, due to which the random walks on it can be 
mapped to the walks on TSF networks 0. Therefore, 
our general conclusion should not be applied to the de- 
terministic LSF networks with high symmetry. 




FIG. 5: (a) Density decay in the DMS networks with m — 1 
and m = 4. The degree distribution exponent is 7 = 2.5. (b) 
Density decay exponent a (circles) and (5 (squares) defined in 
Eq. 11711 for the DMS networks with m = 1 (open symbols) 
and m = 4 (filled symbols). 



III. PAIR-ANNIHILATION PROCESS 

The pair-annihilation process is a diffusion-limited 
reaction-diffusion process. In this process, each node in 
a given network may be occupied by a particle (denoted 
by A) or empty (denoted by 0). The particles perform 
random walks on the network, and annihilate pairwise 
whenever they meet at a same node (A + A — ► 0). 

The pair-annihilation process on SF networks was 
studied numerically by Gallos and Argyrakis They 
found that the particle density decays algebraically 
p(t) ~ t~ a with the 7-dependent decay exponent a = 
ct{l) > L This is contrasted with the d-dimensional pe- 
riodic lattice case where a = a(d) < 1. Namely, the 
particle density decays faster in SF networks. 

Later on, Catanzaro et al. [ill ] developed a mean field 
theory for the pair- annihilation process, which will be 
reviewed briefly hereafter. Let us define pk as the av- 
erage particle density at nodes with degree fc. It is re- 
lated to the total density through the relation p(t) — 
J^fc Pdeg{k)pk{t). In a mean field level, one can write 
down the rate equation for p^ as 

^- = -p k (t) + ^[l-2p k (t)]p(t), (14) 

where (fc) is the mean degree. Multiplying Pd eg (k) and 
summing over fc, one finds that 



dp(t) 
dt 



-2p(t) 



1 

(fc) 



fcP deg (fc) j3fc(f) 



(15) 



Then, Catanzaro et al. made a quasistatic approximation 
neglecting the time derivative in Eq. I|14f> . The approx- 
imation assumes that the particles adjust themselves so 
efficiently that their distribution pu remains close to a 
stationary one at a given value of pit). It leads to the 
relation 



Pk(t) 



k P (t)/(k) 
1 + 2kp{t)/{k) 



Substituting it in Eq. 115(1 and solving the resulting equa- 
tion, Catanzaro et al. obtained that the density decay 
follows a power law in the N — > 00 limit as 



pit) ~t- a (\nty p . 

The decay exponent a is given by 

*)-{; /(7 - 2) :^r 3 



(17) 



(18) 



(16) 



and the exponent for the logarithmic correction is given 
by (3 = 1 for 7 = 3 and f3 = for 7 ^ 3. This result 
is qualitatively consistent with the simulation results of 
Refs. HHEl. 

In the previous works [ToL ITTI] , only the 7-dependent 
scaling behaviors have been studied. However, the study 
on the random walks in Sec. ITTI suggests that the scaling 
behavior may also depend on the global structure of the 
underlying SF network. In this section, we present the 
results of our numerical works on the pair-annihilation 
process on the DMS networks with m = 1 (tree structure) 
and m/1 (looped structure) . Comparing the two differ- 
ent cases, we will show that the global structure does also 
matter for the scaling behavior of the pair-annihilation 
process. 

We have performed the Monte Carlo simulations on the 
DMS networks of size N = 1024000 with various values of 
7 and to. The fully-occupied state is taken as the initial 
configuration. In Fig. 0(a), the numerical data from the 
TSF networks (to = 1) and the LSF networks (to = 4) 
are compared. The data show that the particle density 
decay follows the power law in both cases but with a 
different exponent. 

We estimate the decay exponent by fitting the data to 
the form in Eq. i|17|) . and plot a and f3 as a function of 
7 in Fig. 0(b). At m = 4, our result is qualitatively 
consistent with the previous works; a ~ 1 for 7 > 1, and 
a varies with 7 for 7 < 3. And the logarithmic correction 
is prominent at 7 = 3. However, at m = 1, we obtain a 
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k k 

FIG. 6: pk vs. tc at different time steps t =2° (top),- • -,2 s (bot- 
tom). The data are taken from the DMS network with 
m = 4 (a) and m = 1 (b). The network size is N = 1024000 
and the degree distribution exponent is 7 = 2.5 in both plots. 
The dashed lines have the slope 1. 



However, the particle distribution deviates from 
Eq. (UHJl on the TSF network. Figure ED (b) shows that 
there are three different regimes: pk ~ k for k < ki, 
Pk ~ fc 7 ' for fci < k < k2, and ~ 1/2 for &2 < fc< The 
scaling exponent 77 in the intermediate regime is found 
to be less than 1 and to vary with 7. This feature is 
inconsistent with the quasistatic approximation. 

The quasistatic approximation assumes that particles 
rearrange themselves quickly upon the change in the to- 
tal particle density via diffusion. It requires that the 
diffusion should be a fast process. In the previous sec- 
tion, we have shown that the diffusion is a slow process 
in the LSF networks. This explains why the quasistatic 
approximation is invalid in the TSF network whereas it 
is valid in the LSF network. 

IV. SUMMARY 



completely different result; a ~ 1.0 at all values of 7 and 
the logarithmic correction is present at all values of 7. 
These numerical results show that the global structure 
of the underling SF networks affects the scaling behavior 
of the pair-annihilation process. 

We speculate the origin for the different scaling behav- 
iors. On the LSF networks, the scaling behavior seems 
to be consistent with the analytic mean field result in 
Eq. I|18|) . In the analytic approach one adopts the qua- 
sistatic approximation [ll| leading to the particle distri- 
bution given by Eq. Q16f) . It can be rewritten as 



Pk(t) 



kp{t)/(k) , for k < (k)/p(t) , 
1/2 , for k » (k)/p(t) . 



(19) 



Note that pk oc fc for small k. 

We investigate the particle distribution numerically. In 
Fig- ED (a), we plot Pk(t) against k for the DMS network 
with the looped structure (m = 4). We find that the par- 
ticle distribution pk(t) is fully consistent with Eq. (|19fl . 
This supports that the analytic approach is appropriate 
on the LSF networks. 



In summary, we have investigated the scaling proper- 
ties of the random walk and the pair-annihilation pro- 
cesses on non-deterministic SF networks with the tree 
structure and the looped structure, respectively. In the 
random walk process on TSF networks, we find that the 
relaxation time scales as r ~ N with the network size N 
and that the return probability decay follows the power 
law with the node-dependent exponent. The LSF net- 
work does not display the power-law scalings. In the 
pair-annihilation process, we find that the exponent de- 
scribing the particle density decay is different in TSF and 
LSF networks. Our results show that the global struc- 
ture of the SF network, as well as the degree distribution 
exponent, is the important ingredient in understanding 
the dynamic scaling behaviors. 
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